PROCEEDINGS 

OF SCIENCE 



PoS (LAT2006) 021 
BNL-NT-06/40 



O 
O 

(N 

> 

O 



Lattice QCD at finite density 



CM 
> 



Christian Schmidt* 

Brookhaven National Laboratory, Physics Department, Upton, NY 11955, USA 



E-mail: cschmidt@bnl.gov 



O 



I discuss different approaches to finite density lattice QCD. In particular, I focus on the struc- 
ture of the phase diagram and discuss attempts to determine the location of the critical end-point. 
Recent results on the transiton line as function of the chemical potential (T c (ji q )) are reviewed. 
Along the transition line, hadronic fluctuations have been calculated, which can be used to charac- 
terize properties of the Quark Gluon plasma and eventually can also help to identify the location 
of the critical end-point in the QCD phase diagram on the lattice and in heavy ion experiments. 
Furthermore, I comment on the structure of the phase diagram at large n q . 



X 



XXIVth International Symposium on Lattice Field Theory 
July 23-28, 2006 
Tucson, Arizona, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Lattice QCD at Finite Density 



Christian Schmidt 



T 



-190 
MeV 



quark-gluon 
plasma 




deconfined, 
X-symmetric 



hadron gas 

confined, 

X-SB 



\1 few times nuclear (J, 



matter density 



Figure 1: Sketch of the QCD phase diagram. 



1. Introduction 

Lattice QCD currently is the only quantitative approach to finite temperature QCD based on 
first principle calculation. For a recent review see [[Tj]. At non zero density however, lattice QCD is 
harmed by the sign problem ever since its inception. The Fermion matrix becomes complex and can 
not be interpreted as a probability distribution. Hence straight forward Monte Carlo simulations 
become impossible. For a detailed description of the sign problem in the epsilon regime see [g]. 

During the last few years a lot of progress has been made to circumvent the sign problem 
for small values of \l q jT , where \i q is the quark chemical potential and T the temperature. This 
progress helps to understand the physics relevant for heavy ion collisions and eventually will clarify 
the existence/location of the critical end-point in the QCD phase diagram. In Fig. p] a sketch of the 
QCD phase diagram in the T-jj. plane is shown. Lattice QCD calculations provide more and more 
evidence that the QCD transition at \i q = is not a phase transition in the thermodynamic sense, 
but a smooth crossover. Further evidence was seen recently in [^]. Nevertheless, one can define a 
transition temperature T c by the peak position of the chiral susceptibility. As a function of the quark 
chemical potential the line of transition temperatures {T c {pL q )) is smoothly connected to a critical 
end-point in the (r,/i 9 )-diagram. For larger chemical potentials the QCD transition is expected to 
be first order. At high densities, several color superconducting phases are expected. 

The rest of the article is organized as follows: in Sec. ^|l will briefly recall different methods 
which have been used so far to calculate thermodynamic observables at non zero chemical potential. 
In Sec. |3|l will summarize current knowledge about the /^-dependence of the critical temperature. 
I will continue with reviewing results on quark number fluctuations along the transition line (Sec. |j) 
and the critical point (Sec. ||). Finally I will discuss the physics beyond the critical point in Sec. ^. 

2. Methods to extract information on the chemical potential dependence 

2.1 Reweighting from the {fx = 0)-ensemble 

On the lattice one has to choose several parameters to characterize a thermodynamic system. 
In addition to the number of lattice points in spacial and temporal directions, N s ,N t respectively, we 
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have to choose quark masses m q , the coupling /3 = 6/g 2 and the chemical potential \i q . Together 
these parameters define the lattice spacing a and thus also the temperate T = 1/ aN t and volume 
V = (ctN s ) 3 of the simulated system. A thermodynamic observable is calculated on the lattice as 

<°Wm, = Z(VX) J m ° { detM ( U > m v^ q )} Nf/4 exp{-/3S G (£/)} . (2.1) 

Here Nf is the number of dynamical fermions. The notation is written down for staggered fermions 
as the additional factor of 1/4 in the power of the fermion determinant indicates. See [|J] for a 
discussion of the "4th root trick" needed to calculate the staggered fermion determinant. 

In principle it is possible to calculate the expectation value of the observable at the parameter 
set p = {p,m q ,jj, q }, from an ensemble generated at po = {j6o,»i 9 0,^,0 }• We have the identity 

(0) p = (OR(U;p,po)) po /(R(U;p,p )) po , (2.2) 

where we define the reweighting factor R as 

R(U;p,p Q ) = [detM(U;p)/detM(U; Po )] N f /4 exp{-(/3 - A>)S G } . (2.3) 

The reweighting method as a tool to perform extrapolation and interpolation in the gauge coupling 
j3 goes back to [Bj. For reweighting in the chemical potential it was first used by the Glasgow 
group [01. However, since the overlap between the generated ensemble at /x^o = and the target 
ensemble at \l q > exponentially decreases with increasing \i q , the method was successful only 
after it was generalized to a multi -parameter approach For N t = 4 lattices it was found that 
reweighting along the transition line T c (ix q ) works quite well up to ajj, q <03 or equivalently for 
V q /T<1.2. 

In general the reweighting approach requires the evaluation of the fermion determinant on ev- 
ery configuration. As this is computationally demanding, one may consider to expand the reweight- 



ing factor given in Eq. (2.3) in terms of the chemical potential [pp. In this case the reweighting 



procedure is, however, only correct up to a certain order in H q /T. 

2.2 The Taylor expansion method 

It is conceptually very simple to calculate the expansion coefficients of any observable O 



(Eq. (2.1)) in a Taylor series around \i q = 0: 

O{fl)=c + c l fl + ^c 2 fl 2 + --- . (2.4) 

Since on the lattice all quantities are given in units of the lattice spacing (a), the expansion param- 
eter is ji = ajl q = Nt~ l (fl q /T). This idea goes back to the first calculation of the quark number 



susceptibility [Bfl. The response of meson masses [10] as well as the pressure and further bulk ther 



modynamic quantities [11, 12, 13, 0] have been studied by this method. The first two nontrivial 



coefficients in Eq. ( |2.4| ) are given by 



.'dO\ / dlndetM\ /0 rx 

,'d 2 0\ „ /d<9dlndetM\ / <9 2 lndetM\ / \ /d 2 lndetM\ 
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Besides derivatives of the observable itself, the calculation of derivatives of IndetM with respect 
to jX is required. The derivatives have to be taken at jXo = 0. Note that due to a symmetry of the 
partition function (Z(jU 9 ) = Z(—pL q )) all odd coefficients in Eq. (2.4) vanish identically. For the 
same reason we have (d lndetM/d/t) = at /t = 0. We explicitly use this property in Eq. ( |2.5[ ) to 
derive the expansion coefficients. 

The advantages of this method are that expectations values only have to be evaluated at jX = 0, 
i.e. calculations are not directly affected by the sign problem. Furthermore, all derivatives of the 
fermion determinant can be expressed in terms of traces by using the identity IndetM = TrlnM. 
This enables the stochastic calculation of the expansion coefficients by the random noise method, 
which is much faster than a direct evaluation of the determinant. Moreover, the continuum and 
infinite volume extrapolations are well defined on a coefficient by coefficient basis. 

On the other hand it is a priori not clear for how large n/T the method works and how large 
the truncation errors are. Furthermore one is strictly limited by phase transitions, since phase 
transitions are connected with discontinuities or divergences. An estimation of the convergence 
radius of the series gives a lower bound on the applicability range and thus also a lower bound to 
the phase transition line in the (T,fi) plane (see the discussion in Sec. |5]). 



2.3 Analytic continuation 

At imaginary chemical potentials, the fermion determinant is real and positive, thus simula- 
tions by standard Monte Carlo techniques are possible. Results on the imaginary jXj axis can be 
analytically continued to the real /t« axis. It is especially easy to convert a Taylor series in jXj, 
expanded around jX = 0, into a Taylor series in JXr. Since the series has only even powers of jX, due 
to the the symmetry Z(jX) = Z(—jX), one only has to switch the sign of every second coefficient 
(C2 — > — C2,C(, —> — C6,...). There is however another symmetry of the partition function which 
limits the analytic continuation. Due to the periodicity [ jl"5| ] Z(hr,Hi) = Z(hr,Hi + 2nT/3) simu- 
lations with Hi > will only give access to the physical region (j,r<kT/3. This method was used 



to map out the phase transition line Q16|-[^q]. One should note, that for this method neither an 
evaluation of the determinant, nor any of its derivatives is required. In order to determine the func- 
tional dependence of an observable on /t# one needs many different simulation points for several 
values of jXj, to perform an analytic continuation using a certain Ansatz. 

A demonstration of the imaginary chemical potential method is given in Fig. The method 
was applied to the case of 2-color QCD [|l9|]. Here it is not only possible to calculate observables at 
jX 2 < 0, but also at jX 2 > 0, since 2-color QCD does not suffer from a sign problem. Thus one can 
explicitly check how far the extrapolation in jX 2 is valid and which type of Ansatz is especially well 
suited for the extrapolation. In general, extrapolations with rational functions seem to be better 



than extrapolations with polynomials [|19|, 20]. 



2.4 The canonical approach 

The canonical partition function (Zp) can be constructed by introducing a c) -function into the 
grand canonical partition function which fixes the net number of quarks present in the system. 
Using an integral representation of the 8-f unction one recognizes the integration parameter as an 
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Figure 2: The chiral condensate as a function of p, 2 , Figure 3: Relation between the Baryon number and 
in SU(2) gauge theory [19]. The solid lines are fits the chemical potential from the saddle point approx- 
with different Ansatze to the data for jX 2 < 0. imation, for different temperatures [27]. 



imaginary chemical potential. One finds 

Z C (T,Q = 35) = ^ f_A (|Q exp [-BB^)zGc{T,m) ■ (2-6) 

Thus, the canonical partition functions are the coefficients of the Fourier expansion of the grand 
canonical partition function (Zcg) in imaginary jj.. Here we have used the 2nT /3 periodicity of 
the grand canonical partition function which as a consequence leads to the fact that the canonical 
partition functions vanish for non-integer values of the baryon number B = 3Q. 

The Fourier-coefficients can be computed exactly [^(J. As for the reweighting method (|Q|), 
the evaluation of the fermion determinant is required on every configuration. In fact the same 
method can be used, which is the calculation of all 6N^ eigenvalues of the so-called "reduced 
matrix" [^]. 

After having calculated the canonical partition functions, a relation between the chemical po- 
tential and the baryon number is needed, in order to explore the phase diagram in the (T,/j,) -plane. 
Such a relation can be obtained by using the saddle point approximation of the fugacity expansion 
(which is exact in the thermodynamic limit). One finds 

where p =B/V is the baryon number density and f(T,p) = — ylogZc(T,p) is the Helmhotz free 
energy density. For several different temperatures the saddle point approximation is shown in Fig. || 
[p7|]. Due to the computational costs the calculation have been performed on 6 3 x 4 lattices, with 
Nf = 4 flavors of staggered fermions. As can be seen from the "S"-shape of the curves, one can 
have more than one solution when solving for the baryon density p, at given /I and T. This reflects 
the nature of the transition in the four-flavor theory, which is of first order. By using a Maxwell 
construction, one is able to calculate the two densities p\{T) < PiiJ), giving the lower and upper 
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boundary of the coexistence area, as well as the critical chemical potential Hg r "(T), as functions of 
the temperature. 

2.5 The density of states method 

An alternative to the importance sampling technique used in most Monte Carlo simulations 
is the density of states method. Here one reorders the path integral representation of the partition 
function in the following way: first expectation values with a constrained parameter will be cal- 
culated. I.e. one selected parameter (<p) is fixed. Expectation values according to the usual grand 
canonical partition function (Zqc) can then be recovered by the integral 

<0>=Jd<!> (Of(U))^p((t>) I jd$ (f(u))*P(<i>) ( 2 -8) 

where the density of states (p) is given by the constrained partition function: 

p(x)=Zf(x) = J ®Ug(U)Stf-x) . (2.9) 

Here ( )^ denotes the expectation value with respect to the constrained partition function. In ad- 
dition, the product of the weight functions f,g has to equal the correct measure of Zqc'- fg = 
detM • exp{— Sg}- The idea of reordering the partition functions is rather old and was used suc- 



cessfully for gauge theories [21] and QED with dynamical fermions [22]. For QCD the parameter 



is usually chosen to be the plaquette [23]: <j) = P. In [24], however, the DOS method was con- 
structed for the complex phase (0 = 9). Within the random matrix model, the authors of [25 ] used 
the quark number density (0 = n q ). 

The advantages of this additional integration becomes clear, when choosing (j) =P and g(U) = 
1. In this case p(0) is independent of all simulation parameters. The observable can be calculated 
as a function of all values of the lattice coupling p\ If one has stored all eigenvalues of the fermion 
matrix for all configurations, the observable can also be calculated as a function of quark mass (m) 
and number of flavors^] (Nf). 

Note that this method does not solve the sign problem. It is, however supposed to solve the 
overlap problem. Moreover, it is also possible to combine the DOS method, with the reweighting 



method 2.1, by reweighting the constrained expectation values in the case of g(U) ^ 1. For large 



reweighting distances an overlap problem is then introduced once again. 

3. The transition line 

Using any of the methods presented above, the calculation of the transition line T c (jj, q ) is 
possible. This has been done for many systems, which differ in the number of quark flavors, 
quark masses, physical volume and lattice spacing. This makes a comparison of different methods 
difficult. The case of Nf ■ = 4, m/T = 0.2 and N t = 4, however, has been studied extensively with 
almost all methods. A comparison is given in Fig. |] [27]. As one can see, the agreement between 



different methods is very good up to a\x q as 0.3 or equivalently \X q jT 1.2. For larger chemical 



potential the two results from the reweighting method (Sec. 2. 1 ), indicated as green and blue points, 
seem to stay above the other results. The reason could be the lack of overlap between the simulated 
and the reweighted ensemble. 
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Figure 4: Comparison of the transition line T c (jj.) obtained 
from different methods. Plot from [27]. 
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Figure 5: The mass dependence of the 
transition line for Nf = 3 [29]. 



Other results shown in Fig. |] come from the imaginary chemical potential approach (Sec. 2.3 ), 
solid line, from a generalized imaginary chemical potential approach, black dots, and from the 
canonical approach (Sec. 2.4), red points. These results seem to be in agreement even for somewhat 
larger chemical potentials. However, the transition line from the canonical approach bends down at 
/I /T ~ 1.5. Strong coupling calculations at j3 = [28] show that this is indeed necessary in order 
to match with the correct strong coupling limit. 

To discuss the transition line a bit more quantitatively, one can expand T c (pL q ) in terms of the 
chemical potential 




(3.1) 



In general the first non trivial coefficient ?2 will depend on the number of flavors and the quark 
masses as indicated above. Of course they will also be sensitive to finite volume and cut-off effects. 
One can, however, hope that for large physical volumes and small lattice spacings, i.e. N s /N t > 4 
and N t > 4, those effects are small. A detailed comparison of ?2 is given in Tab. [[} In general, the 
curvature of the transition line becomes steeper for increasing number of flavors and decreasing 
quark masses. Two of the 3-flavor results which have been obtained with the p4-improved action 
[ |29[ ] are shown in Fig. ||. In some cases it has also been possible to estimate the sub-leading 
coefficient ?4, which has been found to be very small or even negative. If an extrapolation with a 
Pade Ansatz is performed, the transition line tends to be steeper for high }X q [ |20| ] compared to the 
truncated Taylor series and shows faster convergence. 

We also note, that the calculation of T c (pL q ) has two parts. The first part involves the calculation 
of j3 c (/l), the second one is the calculation of the lattice /3-function (adfi/da). Some of the results 
listed in Tab. [j] have been obtained by using the perturbative two-loop j8 -function, which has the 
tendency to underestimate the curvature of the critical line. 
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Nr 

J 


am 


N s 


h 


Action 


B -Function 


Method 


Reference 


2 


0.1 


16 


0.69(35) 


p4 


non-pert. 


Taylor+Rew. 


[8] 




0.032 


6,8 


0.500(54) 


stag. 


2-loop pert. 


Imag. 


[17] 


3 


0.1 


16 


0.247(59) 


p4 


non-pert. 


Taylor+Rew. 


[29] 




0.026 


8,12,16 


0.667(6) 


stag. 


2-loop pert. 


Imag. 


[39] 




0.005 


16 


1.13(45) 


p4 


non-pert. 


Taylor+Rew. 


[29] 


4 


0.05 


16 


1.86(2) 


stag. 


2-loop pert. 


Imag. 


[18] 


2+1 


0.0092,0.25 


6-12 


0.284(9) 


stag. 


non-pert. 


Rew. 


[34] 



Table 1: Comparison of the first nontrivial coefficient ti in the Taylor expansion of the transition line. All 
results have been obtained with N t = 4. 



4. Hadronic fluctuations 

Following the transition line into the non-zero chemical potential plane, quark number fluc- 
tuations Xq belong to the most important observables. They will diverge at the critical end-point 
and thus provide an excellent signal for the existence and its location on the lattice and eventually 
may be detectable in heavy ion experiments. Hadronic fluctuations can be computed from Taylor 
expansion coefficients of the pressure with respect to the quark chemical potential: 



^ = i>(r)(^)" with Cn (T) 



1 N? dlnZ 



T 4 to KtJ n\Nfd(fiN t )" 



(4.1) 



Due to the particle anti-particle symmetry {pL q «-> — ji, q ) all odd coefficients vanish. Thus the first 
three non-zero coefficients are c%, c/\, and C(,. They have been calculated in the case of two flavors 
of p4-improved staggered fermions, with m q /T = 0.4 [ jl3"| ] and are shown in Fig. ^. Note that in 
the Taylor expansion of the pressure the up and down quark chemical potentials have been chosen 
to be equal. Having calculated the coefficients c n {T) one can construct the quark number density 
and quark number fluctuations 

S=E^(n(^)"" 1 ; ^ 2 =in(n-l)c n (T)mY- 2 . (4.2) 

1 n=2 K 1 ' 1 n=2 K 1 ' 

Also shown in Fig. |6| are the quark number fluctuations for various values of the chemical potential. 
It is interesting to see that at \i q = 0, the fluctuations % q show a rapid but monotonic increase at the 
transition temperature, whereas a cusp is developing at T c (pL q ) for \i q > 0. This is a clear sign for 
approaching the critical end-point. 

In the case of 2-flavor QCD, the quark number susceptibility is directly proportional to the 
baryon number fluctuation. Having two light quarks and one heavier strange quark, the situation is 
not that simple anymore. To match the situations realized in heavy ion collisions, one still wants to 
expand the pressure in terms of \i q = pL u = jiij, keeping the strange quark chemical potential zero 
(Us = 0). However, in order to analyze fluctuations of conserved quantum numbers, it appears to 
be more appropriate to perform a basis change going from the space of quark number fluctuations 

^ = M-(n a )(n p )o,ll^ , (4.3) 



8 



Lattice QCD at Finite Density 



Christian Schmidt 




where a,j8 6 {u,d,s} to the space of hadronic fluctuations, indicated by a,fi € {It,,Y,B} with Ij, 
being the third component of the isospin, Y being the hypercharge and B the baryon number. In 
Fig. [7| the diagonal susceptibilities ,/ 3 , Xyj and %b,b as well as the off-diagonal susceptibility %y,b 
are shown [30]. They have been measured by the MILC Collaboration with 2+1 flavor of Asqtad 
fermions on a 12 3 x 6 lattice. The light quark mass is m q = 0.2m s where m s is the physical strange 
quark mass. The curves have been normalized such that the continuum Stefan-Boltzmann value 
is 0.5 for all of them. Qualitatively they show the same behavior as the diagonal quark number 
fluctuations, only the off-diagonal susceptibility %y,b shows a cusp already at jj, = 0. However, up 
to a minus sign also the off-diagonal quark susceptibility % u j shows a cusp at /I = 0. 

Furthermore, also fluctuations of other conserved quantities such as electric charge Q have 
been computed. It is especially interesting to analyze the correlations between different conserved 
charges. Also shown in Fig. ^ are the correlations 

_ (XS) - (X) (S) 

Cx > s ~ (si)-(xf ' (4 " 4) 

where S is the strangeness and the operator X is either the electric charge (X = Q) or the baryon 
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Figure 7: Different hadronic fluctuations [30] and the correlation between strangeness and baryon number, 
electric charge fluctuations respectively [31]. 



number (X = B) [^TJ]. Such calculations have been performed on an 16 3 x 4 lattice, using a standard 
staggered action. The two dynamical light quark masses yield m n /m p = 0.3. The strange quark 
has been treated in the quenched approximation. The results show that above T c strangeness and 
electric charge or baryon number fluctuate independently. This is consistent with the quasi-particle 
picture of the Quark-Gluon-Plasma (QGP). However, below T c , in the hadronic phase fluctuations 
are clearly correlated. These correlations are thus directly related to the degrees of freedom in the 
QGP. These are also clearly visible in the calculation of higher order cumulants of fluctuations [[32|]. 



5. The critical end-point 



Locating the critical point is one of the most challenging tasks for lattice QCD at finite chem- 



ical potential. The first attempt to locate the critical point used the reweighting method [|33[]. For 
this calculation, 2+ 1 flavor of standard staggered fermions have been used at a pion mass of about 
300 MeV and a kaon mass of about 500 MeV. Lattice sizes have, however, been rather small (4 3 x 4 
- 8 3 x 4). A critical chemical potential of \i c ^" = 725(35) was found. A second calculation [Ql, 
using again the reweighting method, with physical masses {m n = 150 MeV, = 500 MeV) and 
somewhat larger volume (6 3 x 4 - 12 3 x 4), let to pLg U = 360(4) MeV. 

When using the reweighting method for locating the critical point, the minima of the normal- 
ized partition function in the complex /3 -plane (Lee-Yang zeros) have to be determined 



Z(/3Re,j6i m)i U) 



Z(i8 R e,0,0) 



6rj3W,W J 3 AS Ce (e g (W / /4)(lndetM(^)-lndetM(0))\ 



(5.1) 



/ (&e,0,0) 

1 , this can be done with high accuracy, which can be 



In SU (3) gauge theory, where we have e' e 
seen in Fig. M [35]. One can even identify a second Lee- Yang zero. In order to locate the critical 



point, one has to take the infinite volume limit and monitor the approach of the Lee- Yang zeros 
on the real axis. When the first Lee- Yang zero touches the real axis in the infinite volume limit, 
a critical point has been reached. In Fig. || this is shown for full QCD with physical masses [Q]. 
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Figure 8: Lee-Yang Zeroes in the complex j3 plane, in the case of SU (3) gauge theory [35] (left) and the 
distance of the smallest Lee- Yang zero from the real axis as function of the chemical potential, in the case 
of full QCD [34] (right). 



In QCD with non-zero chemical potential the analysis of Lee- Yang zeros is, however, subtle [J35]]. 
For large volumes and chemical potentials the phase factor of the determinant e' e will force the 
Lee- Yang zero onto the real axis, which might lead to an underestimation of the critical point. 

Another difficulty with the reweighting method at finite chemical potential has been pointed 
out in [36]. It was noted, that taking the fourth (or square) root of the determinant (which is 
necessary in order to simulate 2 or 1-flavor QCD with staggered fermions; see also [Q]) could lead 
to phase ambiguities. This problem becomes acute when \i q > m n /2. 

All of the a above mentioned limitations are, however, irrelevant for the location of the critical 
point with the reweighting method if the critical point is located at small values of \i q . 

Using the Taylor expansion coefficients of the pressure, it is also possible to estimate the 
location of the critical point. The convergence radius of the expansion is limited by the nearest 
singularity in the complex chemical potential plane. For each fixed temperature, the radius of 
convergence is given by 



lim p n 




Cn+2 



(5.2) 



Moreover, the sign of the coefficients c n gives information about the location of the singularity 
in the complex plane. If all coefficients are positive, the singularity is located on the real axis of 
the complex chemical potential plane. If the sign is strictly alternating, the singularity lies on the 
imaginary axis. For a detailed discussion see [37]. 

Having only a limited number of expansion coefficients, one can only estimate p. The hope 
is that the convergence of the p„ will be fast. Indeed, a clustering of the p„ is seen in the phase 
diagram, as shown in Fig. |9] [13]. This calculation, which has been performed with 2 flavors of p4 
improved fermions and m n jm p = 0.7, suggests a critical chemical potential of jig" pa 500 MeV. 



11 



Lattice QCD at Finite Density 



Christian Schmidt 



2.5 



1.5 



0.5 



77T A 



SB(p ) 



'SB(p.) 



e ■ — 

-e — 




□ P„ 
o P 2 
P 4 



i I i I , L_ 

2 4 6 



r(n ) 



~q 



10 12 



0.1 



c 4 /c 2 



n f =2+1, m q =0.1m s 
n f =2, m q /T=0.4 



cont. 



T/T n 





9 


® 


N t =4 







0.8 1 1.2 1.4 1.6 1.8 2 



Figure 9: Estimates of the radius of convergence in the (r,|i ? )-plane (left), the ratio C4/C2 of the expansion 
coefficients (right). 



All calculated p„ are, however, consistent (within statistical error) with the resonance gas model in 
the Boltzmann approximation, where the radius of convergence is infinity. 

The authors of p8[ ] have estimated the critical chemical potential from a Taylor expansion 
of the quark number susceptibility and find \i# lt 180 MeV. Two flavors of standard staggered 
fermions have been used on lattices up to 24 3 x 4 and quark mass corresponding to m n /m p = 0.3. 
The difference between the two estimates [13, 38] of the critical point is large. We note that the 
second estimate comes from the expansion coefficients of % q . As can be seen from Eq. this 
will result in a smaller p„ for each fixed n. The limit limbec p n is of course the same. For finite n, 
however, the estimate of \i c ^" « 1 80 MeV would correspond to \i^ lt « 240 MeV, when estimating 
the p„ with coefficients of the same order from the expansion of the pressure. Nonetheless, the 
difference between the two estimates is still striking. The origin could be the difference in mass. 
However, preliminary results from the RBC-Bielefeld Collaboration, also shown in Fig. ^, do not 
indicate a strong mass dependence in C4/C2 = l/pf- 

The critical point can also be studied directly at \i q = 0. This can be done by tuning the quark 
masses carefully to a value where the critical chemical potential is \i# lt = 0. In the quark mass 
plane of two degenerate light quark and one strange quark, (m u ^ , m s )-plane, a line exists on which 
this condition is fulfilled. Starting from this line, one can define a surface of critical points in the 
3-dimensional space of (m u ^,m s ,^i q ). On one site of the surface, the order of the QCD transition is 
first order (for smaller masses) on the other side the transition is only a smooth crossover. The line 
of critical points at \i q = has been computed for standard staggered fermions [39] as shown in 
Fig. 10. For locating the critical points, the fourth order Binder cumulants of the chiral condensate 
have been calculated. Since the probability distribution of the order parameter is universal, also 
the value of its forth cumulant is a renormalization group invariant quantity which furthermore is 
volume independent at the critical point. From a calculation of cumulants at imaginary chemical 
potential the region of first order phase transitions was found to shrink as sketched in Fig. 1C. 
This calculation has been performed with standard staggered fermions on an rather course lattices 
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E 




Figure 10: Line of critical points in the quark mass plane for fig = and i[iB = 2 AT (left) and sketch of the 
surface of critical points as found in [39] (right). 



(N t = 4). If this result gets confirmed in the continuum limit it would put doubts on the existence 
of a critical point in QCD with physical masses. It is interesting to mention, that a shrinking region 
of first order transitions has also been found in the case if isospin chemical potential pT] ] . 

6. Beyond the critical point 

Even more challenging than locating the critical point, is the study of the physics at high den- 
sities and low temperatures. One attempt to do so is a calculation using the density of states method 
[p0[l. Using four flavors of standard staggered fermions (i.e. taking the root of the determinant is 
no necessary), several simulation points in the (j8,jLt) plane have been chosen to generate phase 



quenched configurations by employing the method proposed in [|42|]. The lattice size has been 
6 3 x 4 and 6 4 . The quark mass was chosen to be m/T = 0.2. The generation has been done with 



constrained plaquettes. In oder to do so, the 5-function in Eq. 2.9 has been replaced by a sharply 
peaked Gaussian potential, which in practice means that the force term in the HMD-R algorithm 
had to be modified. In the notation used in Sec. 2.5 this would mean (j) = P, g = |detM|e _ ^ Sc and 



f = e . For each simulation point, several runs have been performed with about 20 different values 
of the plaquette. By calculating the eigenvalues of the reduced matrix the phase of the determinant 
was calculated for each of those runs. By numerically calculating the integrals 

(P) = j dxxp(x){cos(d)) x (P 2 ) = J dxx 2 p(x) (cos(d)) x , (6.1) 

we recover the grand canonical expectation value of the plaquette and its square. Here p(x) is 
the density of states (Eq. |2.9| ), which has been measured by the integral method, usually used to 
calculate the pressure. The susceptibility of the plaquette is then given by the usual expression 
%p = (P 2 ) — (P) 2 . From the peak position of the plaquette susceptibility the phase diagram was 
calculated as shown in Fig. 1 1 . The scale was set by the Sommer radius ro, measured on a 10 3 x 20 
lattice. We find a triple point, where three different phases seem to coexist. The phases show 
different plaquette expectation values. The triple point is located around /i^ n as 300 MeV, however 
its temperature (r tri ) decreases from T tri « 148MeV on the 6 3 x 4 lattice to T bi sa 137MeV on the 
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Figure 11: Phase diagram from the density of state method [39]. 



6 4 lattice. This shift reflects the relatively large cut-off effects one faces with standard staggered 
fermions and temporal extents of 4 and 6. 



Also shown in Figure |ll] are points from simulations with quark mass m/T = 1.2. The phase 



boundary turned out to be — within our statistical uncertainties — independent of the mass. 
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